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Resumen 


Campos-Aranda, D. F. (mayo-junio, 2015). Búsqueda del 
cambio climático en la temperatura máxima de mayo en 16 
estaciones climatológicas del estado de Zacatecas, México. 
Tecnología y Ciencias del Agua, 6(3), 143-160. 


De manera general y aproximada, el planeta es 0.8 *C 
más caliente que en el siglo XVIII antes de la Revolución 
Industrial y podrá ser 2 *C el aumento en el lapso de una 
generación. En tales circunstancias, resulta trascendental 
cuantificar la magnitud e intensidad del cambio en cada 
una de las variables climáticas básicas. La identificación 
del cambio climático (CC) se debe apoyar en las técnicas 
estadísticas, para buscar en los registros largos tendencias 
que demuestren su presencia de manera significativa. En 
este trabajo se buscó el CC en la temperatura del estado 
de Zacatecas, México, para lo cual se procesaron 16 
estaciones climatológicas con más de 50 años de registro de 
temperatura máxima del mes de mayo. Se aplicaron pruebas 
estadísticas generales y específicas para buscar componentes 
determinísticas. Tres de estas pruebas, las relacionadas con 
la detección de la tendencia y su verificación, se describen 
con detalle en apéndices. Se encontró que cuatro de las ocho 
estaciones climatológicas de la Región Hidrológica 12 Parcial 
(Río Santiago) presentan tendencia ascendente significativa, 
cuyo valor medio del incremento entre la media reciente y la 
de largo periodo es de 2.1 *C. 


Palabras clave: homogeneidad, persistencia, tendencia, 
cambio en la media, pre-blanqueado, tendencia lineal. 


Introducción 
Generalidades 
Los cambios naturales del clima han ocurrido 


durante toda la historia de la Tierra, lo preocu- 
pante ahora es el calentamiento global originado 


Abstract 
Campos-Aranda, D. E. (May-June, 2015). Search for Climate 
Change in Maximum Temperatures in May based on 16 Weather 
Stations in State of Zacatecas, Mexico. Water Technology and 
Sciences (in Spanish), 6(3), 143-160. 


Overall, the planet is roughly 0.8%C warmer than in the 18th 
century, before the industrial revolution, and could increase to 
2*C warmer by the next generation. Given these circumstances, 
it is crucial to quantify the magnitude and intensity of change 
in each of the basic climatic variables. Climate change should be 
determined based on statistical techniques that search long registries 
for trends that significantly demonstrate its presence. This work 
seeks to identify climate change based on temperature in the state of 
Zacatecas, Mexico. To this end, information from 16 weather stations 
was processed, with over 50 years of high temperatures recorded 
for the month of May. General and specific statistical tests were 
applied to search for deterministic components. Three of these tests, 
the ones related to detecting and verifying trends, are described in 
detail in the appendices. Four of the eight weather stations in Partial 
Hydrological Region 12 (Santiago River) presented significant 
increasing trends, with a mean value of the increases between the 
recent mean and the long-term mean of 2.1*C 


Keywords: Homogeneity, persistence, trend, change in mean, 
prewhitening, linear trend. 
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por las actividades humanas, el cual se ha de- 
sarrollado en un tiempo breve. De acuerdo con 
el Cuarto Informe del Panel Intergubernamental 
sobre Cambio Climático (IPCC, 2007), si las 
concentraciones de gases de efecto invernadero 
se mantuvieran constantes con la magnitud que 
tenían en el año 2000, el incremento en la tem- 
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peratura en la década de 2090-2099 relativa a 
la del periodo 1980-1999, será de 0.6 *C, con un 
intervalo probable de 0.3 a 0.9 *C. Sin embargo, 
para el escenario más bajo o favorable de con- 
centraciones (B1), el incremento citado será de 
1.8 *C, con un rango de 1.1 a 2.9 *C y para el 
escenario más alto o crítico (A1FD), los valores 
mencionados serán 4.0 *C, variando de 2.4 a 6.4 
*C. Este calentamiento global es considerable y 
seguramente intensificará el ciclo hidrológico, lo 
cual resultará en un incremento de la frecuen- 
cia y magnitud de las crecientes y las sequías, 
asociado lógicamente con una modificación de 
los patrones de ocurrencia de la temperatura y 
la precipitación (Cavadias, 1995; IPCC, 2007). 

Actualmente, los procedimientos de pla- 
neación, diseño y operación de los aprovecha- 
mientos hidráulicos se basan en la consideración 
de que los procesos hidrológicos cambian de 
forma muy lenta, de manera que pueden ser 
considerados estacionarios. Los únicos efectos no 
estacionarios que han sido tomados en cuenta 
son la evolución de los métodos de medición 
y los cambios en las condiciones físicas de las 
cuencas, originadas por la deforestación, prácti- 
cas agrícolas, urbanización, etcétera (Yevjevich, 
1995). 

Aunque los cambios en las cuencas y en las 
demandas pueden tener un mayor impacto 
en los proyectos de aprovechamiento de los 
recursos hidráulicos que los debidos al cambio 
climático, los primeros son eminentemente loca- 
les y los segundos regionales o nacionales, con 
consecuencias graves durante la vida útil de la 
obra. Por lo anterior, resulta necesario eviden- 
ciar el cambio climático y formular estrategias 
que se anticipen a sus consecuencias (Martínez- 
Austria €: Aguilar-Chávez, 2008). 


Estudios de búsqueda de tendencias climáticas 
en la temperatura 


Arora, Goel y Singh (2005) procesaron los 
registros de temperatura media, máxima y 
mínima de 125 estaciones climatológicas de la 
India, según cinco series de datos, la anual y en 
cuatro épocas: invierno (de diciembre a febrero), 


pre-monzón (de marzo a mayo), monzón (de 
junio a agosto) y pos-monzón (de septiembre 
a noviembre). Usaron el test no paramétrico de 
Mamn-Kendall para detectar la tendencia y en- 
contraron que ésta es ascendente en el lapso de 
1941 a 1999, con valores promedio de 0.42, 0.92 
y 0.09 *C por cada 100 años, en la temperatura 
anual, máxima y mínima, respectivamente. En 
el ámbito regional, en el sureste y oeste de la 
India se encontró tendencia ascendente y en el 
norte, descendente. Los porcentajes de tenden- 
cias ascendentes o descendentes significativas 
en cada época y para cada temperatura varían; 
por ejemplo, para la temperatura máxima en la 
época de pre-monzón se tienen 18 estaciones 
con tendencia descendente y 42 con ascendente; 
entonces existe un 48% de tendencias signifi- 
cativas. 

Singh, Kumar, Thomas y Arora (2008) anali- 
zaron los registros de 90 a 100 años de tempera- 
tura media, máxima y mínima de 43 estaciones 
climatológicas de las regiones noroeste y central 
de la India, las cuales abarcan nueve cuencas de 
grandes ríos. Usaron el test no paramétrico de 
Mann-Kendall para detectar la tendencia, usan- 
do series pre-blanqueadas si existió persistencia 
(ver Apéndice 3). En siete cuencas encontraron 
tendencia ascendente en la temperatura media 
con valores que variaron de 0.40 a 0.60 *C por 
cada 100 años y en dos cuencas la tendencia fue 
descendente con magnitudes de -0.15 a -0.44 
*C por cada 100 años. En el área de estudio, la 
temperatura máxima mostró tendencia creciente 
con una magnitud media de 0.81 *C por cada 
100 años; este valor resultó mayor que el obser- 
vado para la temperatura mínima. 

Martínez, Maleski y Miller (2012) buscaron 
tendencias en las temperaturas medias, máxi- 
mas y mínimas del estado de Florida, Estados 
Unidos, analizándolas por meses, épocas y sus 
valores anuales. Con el test de Mann-Kendall 
procesaron 22 estaciones climatológicas y 
emplearon dos lapsos de registro, 1985-2009 
y 1970-2009; encontraron para la temperatura 
máxima en ambos intervalos de análisis y en 
todas las duraciones empleadas, que unas esta- 
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ciones muestran tendencias ascendentes y otras 
descendentes, con cierta coherencia espacial. 

Sonali y Nagesh-Kumar (2013) estudiaron 
las series integradas de temperatura máxima y 
mínima en cada una de las siete regiones ho- 
mogéneas que definieron para la India, usaron 
valores anuales, mensuales y en cuatro épocas: 
invierno (enero y febrero), pre-monzón (de mar- 
zo a mayo), monzón (de junio a septiembre) y 
pos-monzón (de octubre a diciembre); así como 
tres intervalos (años) de análisis: 1901-2003, 
1948-2003 y 1970-2003. En la detección de la 
tendencia emplearon seis técnicas no paramé- 
tricas. Encontraron una tendencia ascendente 
consistente en las siete regiones homogéneas 
durante las tres últimas décadas, con una mayor 
pendiente para la temperatura mínima que para 
la máxima. 


Proyecciones climáticas para el estado de 
Zacatecas 


De acuerdo con Montero-Martínez y Pérez- 
López (2008), para el estado de Zacatecas, el 
incremento en la temperatura media de su- 
perficie en el verano (julio-agosto-septiembre) 
en la década de 2050 será de 1.8 a 2.5 *C; en el 
mismo lapso, pero en 2070, de 2.5 a 3.2 *C, y en 
la década de 2090 de 4.0 a 4.7 *C. En específico, 
para la temperatura máxima del verano (junio- 
julio-agosto) en Zacatecas, Montero-Martínez, 
Martínez-Jiménez, Castillo-Pérez y Espinoza- 
Tamarindo (2010) definen los 
siguientes para el periodo de 1961-1990: en la 


intervalos 


zona centro del estado varía de 24 a 27 *C; en su 
región norte fluctúa de 30 a 33 “C, y en las áreas 
sur y media (entre centro y norte) cambia de 
27 a 30 *C. En cuanto a las proyecciones de sus 
anomalías promedio, Montero-Martínez et al. 
(2010) pronostican para el lapso de 2010 a 2039 
un aumento general en el estado de Zacatecas 
de 1.0 a 1.5 *C; para el intervalo de 2040 a 2069, 
el incremento será de 2.5 a 3.0 *C; por último, 
para el periodo de 2070 a 2098, encuentran un 
aumento en las zonas centro y norte de 4.0 a4.5 
*C, y en sus áreas sur y media de 4.5 a5.0 *C. 


Montero-Martínez et al. (2010) también 
indican que la regionalización (downscaling) de 
las proyecciones climáticas generadas con los 
modelos climáticos globales acoplados es una 
parte fundamental de todo estudio orientado a 
cuantificar el probable cambio climático, el cual 
constituye el llamado enfoque dinámico. La otra 
técnica de trabajo en la búsqueda del cambio 
climático se denomina enfoque estadístico; en ella 
se utilizan registros climatológicos históricos a 
los cuales se les aplican diversos algoritmos o 
técnicas estadísticas (Mujumdar € Nagesh- 
Kumar, 2012). Este segundo enfoque tiene dos 
áreas de análisis: el pasado y la proyección hacia 
el futuro. En la acepción del estudio del pasado, 
la variabilidad climática puede generar una apa- 
rente tendencia cuando se analiza un registro 
corto; por ello, Kundzewicz y Robson (2004) 
establecen que en los análisis de detección del 
cambio climático se utilicen registros de cuando 
menos 50 años. 


Objetivo del estudio 


En este trabajo se buscan cambios significativos 
estadísticamente en la temperatura máxima 
mensual del mes más caluroso, el cual resultó 
ser mayo, en las estaciones climatológicas de 
amplio registro (> de 50 años) del estado de 
Zacatecas, México, con el objetivo de ilustrar 
cómo se deben procesar los registros de la 
variable climática temperatura, sea ésta media, 
máxima o mínima, y sean sus series de datos 
anuales, mensuales o por épocas. Se procesaron 
16 registros, cuyo número de datos varió de 50 
a 67 años. Se aplicaron tres pruebas estadísticas 
generales para detectar pérdida de homogenei- 
dad y 11 específicas para buscar componentes 
determinísticas, como persistencia, tendencia 
y/o cambios en la media o en la variabilidad. 
Todos los tests estadísticos relacionados con la 
búsqueda de la tendencia ascendente y otros 
procedimientos asociados con su verificación se 
exponen con detalle en apéndices. Se encontró 
que cuatro estaciones climatológicas de las ocho 
localizadas en la Región Hidrológica 12 Parcial 
(Río Santiago) del estado de Zacatecas presentan 
tendencia ascendente significativa. 
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Materiales y métodos cuadro 1, por región hidrológica y en orden 
decreciente de altitud. En la figura 1 se muestra 
Estaciones climatológicas seleccionadas su ubicación geográfica dentro del estado de 


Zacatecas, México. 


Este trabajo busca detectar los efectos del 


cambio climático en la temperatura en el estado Registros de temperatura máxima procesados 
de Zacatecas, México, y para ello se consideró 

conveniente estudiar la temperatura máxima Con base en los registros mensuales de tem- 
del mes más caliente, pues seguramente en peratura máxima se buscó el mes más cálido 
estos registros serán evidentes tales efectos. tanto por su promedio mensual como por sus 
El registro en Excel de temperatura máxima valores máximo y mínimo; se encontró que en 
(*C) mensual de la Dirección Local Zacatecas la mayoría de las estaciones climatológicas se- 
de la Comisión Nacional del Agua (Conagua) leccionadas, tal mes es mayo; sólo en la estación 
cuenta con 133 estaciones climatológicas, de El Platanito resultó ser junio. También en varias 
las cuales sólo 27 tienen registros de más de estaciones se encontró que junio tiene un valor 
50 años, pues comenzaron a operar en 1963 mínimo más alto que mayo, como en Chalchi- 
O antes. Sin embargo, de tales estaciones, 11 huites, El Sauz, La Florida, San Pedro Piedra 
tienen lapsos de registro ausente y/o varios Gorda y Juchipila. En la estación El Cazadero, 
años incompletos debido a datos mensuales mayo y junio tienen valores iguales del prome- 
faltantes. Por lo anterior, sólo se pudieron dio y sus extremos. Los 16 registros procesados 
obtener 16 estaciones con registros mayores corresponden a la temperatura máxima del 
de 50 años y con pocos años incompletos, mes de mayo, cuyos escasos valores faltantes 
cuyas características generales se citan en el se tomaron igual a su media. 


Cuadro 1. Características generales de las 16 estaciones climatológicas procesadas del estado de Zacatecas, México. 


Núm Nombre as Els Periodo de registro (NA?) 
msnm núm. 
1 Chalchihuites 2260 11 1962-2012 (51) 
22 Presa Santa Rosa 2240 36 1947-2012 (66) 
3 El Sauz 2 090 36 1947-2012 (66) 
4 El Cazadero 1920 36 1959-2012 (54) 
5 Zacatecas 2 485 37 1953-2012 (60) 
6 Calera 2192 37 1958-2012 (55) 
7 Agua Nueva 1932 37 1963-2012 (50) 
e 8 Monte Escobedo 2190 12 1963-2012 (50) 
3 9 Villa García 2120 12 1950-2012 (63) 
= 10 San Pedro Piedra Gorda 2 032 12 1949-2012 (64) 
El 11 La Florida 1870 12 1955-2012 (58) 
E 12 La Villita 1790 12 1958-2012 (55) 
no 13 Excamé 1740 12 1946-2012 (67) 
á 14 Presa El Chique 1620 12 1950-2012 (63) 
= 15 Juchipila 1270 12 1947-2012 (66) 
E 16 El Platanito 990 12 1958-2012 (55) 
> ¡Región hidrológica. 


"Número de años. 
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Figura 1. Localización geográfica de las 16 estaciones climatológicas procesadas del estado de Zacatecas, México. 


Pruebas estadísticas aplicadas 
Generalidades sobre requerimientos estadísticos 


Cuando los registros de lluvia, temperaturas o 
escurrimiento mensuales o anuales son utiliza- 
dos para obtener parámetros hidrológicos, que 
son la base de los diseños de las obras hidráuli- 
cas de aprovechamiento, tales series o muestras 
deben tener ciertas características estadísticas, 
para que dichas estimaciones sean confiables. 
En resumen, se requiere que las series sean 


homogéneas o consistentes, libres de tendencia 
y que se hayan originado por un proceso es- 
tocástico, cuya componente aleatoria siga una 
determinada distribución de probabilidades 
(Adeloye $: Montaseri, 2002). 

La consistencia implica que los datos o mues- 
tra procedan de una sola población, por ello se 
dice que sean homogéneos. La tendencia existe 
cuando los datos muestran una correlación 
(positiva o negativa) apreciable con el tiempo, 
de manera que su media aritmética no es cons- 
tante y por ello la serie deja de ser estacionaria. 
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La tendencia se origina por modificación de 
las condiciones circundantes a la estación de 
medición, como son los cambios en el uso del 
suelo, o bien por el cambio climático. Relacio- 
nada con la tendencia, se tiene la persistencia 
que caracteriza la propensión de los datos de 
una serie al recordar la magnitud de sus valores 
antecedentes o a ser influenciados por ellos. 
Debido a lo anterior, se cuantifica con base en 
el coeficiente de correlación serial que define la 
dependencia entre cada elemento i y el ¡-k; por 
ello k se llama “tiempo de retraso” (Adeloye $ 
Montaseri, 2002; Machiwal é Jha, 2008, 2012). 


Pruebas y test aplicados 


Se aclara que no existe diferencia semántica en- 
tre prueba y test, ambos estadísticos. Existen dos 
tipos de pruebas estadísticas: las paramétricas 
y las no paramétricas; las primeras requieren 
que los datos provengan de una distribución 
normal, son más potentes o eficientes que las 
segundas, las cuales han adquirido populari- 
dad debido a que no presentan tal restricción. 
Machiwal y Jha (2012) han establecido que para 
obtener la misma conclusión con igual nivel de 
significancia (a), las pruebas no paramétricas 
requieren de un 5 a un 35 % de más datos que 
los test paramétricos. 

Debido a lo anterior, de inicio se determinó 
con los test del Cociente de Geary (CG) y de 
Shapiro-Wilk (SW) si en los registros procesados 
de temperatura máxima del mes de mayo se 
puede aceptar que provienen de una distribu- 
ción normal. Los resultados de ambos test se 
tienen en el cuadro 2, en sus columnas 3 y 4, 
indicando, para el segundo, la normalidad con N 
y la no procedencia de una distribución normal 
con NN. 

Las pruebas estadísticas aplicadas se dividen 
en dos grupos: generales y específicas. Las pri- 
meras fueron las tres siguientes: (1) de Helmert 
(PH); (2) de las Secuencias (PS), que detectan 
si la serie O registro es homogéneo debido a que 
no tiene persistencia o tendencia, cambio en la 
media o bien no oscila demasiado; (3) la prueba 


de Von Neumann (PV) es de aleatoriedad contra 
componentes determinísticas no especificadas. 
Los resultados de estas pruebas se exponen en 
las columnas 5 a 7 del cuadro 2, con H para 
registros homogéneos y con NH para los no 
homogéneos. 

Respecto a los test específicos, se aplicaron 
los 11 siguientes: (1) dos de persistencia, el de 
Anderson (TA) y el de Sneyers (TS); (2) tres de 
tendencia, el de Kendall (TK), Spearman (TS,) y 
Mann-Kendall (TM); (3) uno de variabilidad, el 
de Bartlett (TB); (4) uno de cambio en la media, 
el de Cramer (TC), y (5) cuatro que buscan cam- 
bios en la media o pérdida de homogeneidad 
con estadísticos basados en la curva masa resi- 
dual (TB, a TB,). Todas las pruebas se aplicaron 
con un nivel de significancia (a) de 5%. En el 
Apéndice 1 se citan las referencias donde se 
pueden consultar 13 de los test anteriores, pues 
los tres relativos a la tendencia se exponen en 
el Apéndice 2. Los resultados de estos 11 test 
específicos se muestran en las columnas 8 a 18 
del cuadro 2, usando H cuando no detectan 
la componente que buscan y NH cuando la 
encuentran. 


Análisis de los resultados 
Respecto a las pruebas estadísticas 


El test de Shapiro-Wilk detecta que en cinco de 
los registros procesados de temperatura máxi- 
ma del mes de mayo no se puede aceptar que 
proceden de una distribución normal. En tales 
registros se obtuvieron valores de cociente de 
Geary tan alejados de la unidad como 0.853 y 
tan cercanos como 0.960 (ver cuadro 2). 

Con los test generales se encuentra que sólo 
dos registros, El Cazadero y Monte Escobedo, 
son homogéneos; el resto resulta no homogé- 
neo, la mayoría en las tres pruebas. "Todos los 
registros presentan persistencia y sólo en cuatro 
es detectada exclusivamente por la prueba de 
Sneyers, pudiéndose decir que El Cazadero, 
El Platanito, Monte Escobedo y Excamé tienen 
persistencia ligera. 
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Respecto a la tendencia, seis registros la 
presentan, cuatro con las tres pruebas aplicadas 
y expuestas en el Apéndice 2: Monte Escobedo, 
San Pedro Piedra Gorda, La Villita y Juchipila; 
Villa García sin tendencia, según prueba de Ken- 
dall, y Excamé con tendencia de acuerdo con la 
prueba de Spearman. Con excepción de Excamé, 
el resto de registros que mostraron tendencia 
presentan evidencia de cambio en la media se- 
gún las pruebas (TB, a TB,) de Buishand (1982), 
basadas en la cuerva masa residual, por lo cual 
es probable que su tendencia sea significativa. 

Sólo tres registros mostraron falta de homo- 
geneidad con respecto a la variabilidad (Test de 
Bartlett): Chalchihuites, que tiene dos valores 
muy bajos, y Presa El Chique y Juchipila, que 
por el contrario tienen magnitudes muy altas 
registradas (46.5 "C). Sin embargo, tales resul- 
tados no son confiables, pues esos tres registros 
no proceden de poblaciones normales. 


Verificación de la tendencia detectada 


En los seis registros que mostraron tendencia 
con alguno de los test se aplicó la prueba de 
Mann-Kendall con la serie pre-blanqueada y con 
el pre-blanqueado libre de tendencia, cuyos pro- 
cesos operativos se detallan en el Apéndice 3; 
los resultados se han concentrado en el cuadro 
3, observando que sólo el registro de la estación 
Excamé no presenta tendencia en las pruebas 
de Mann-Kendall aplicadas a su serie preblan- 
queada (X') y con pre-blanqueado sin tendencia 
(Y”). También a los seis registros con tendencia 
se les ajustó un modelo de regresión lineal 
para representarla y se buscó si su pendiente 
resultaba estadísticamente diferente de cero con 
una prueba basada en la t de Student, la cual se 
expone en el Apéndice 4; los resultados se citan 
en el cuadro 3. Exclusivamente en el registro de 
la estación Excamé su pendiente es igual a cero. 

Con el objeto de verificar numéricamente 
las pendientes (1m) estimadas con la regresión 
lineal (Apéndice 4), que se sabe son afectadas 
por los valores dispersos, se aplicó el criterio 
(ms) de Sen (1968) expuesto en el Apéndice 5. 
Sus resultados se tienen en la última columna 


del cuadro 3. Se observa que los valores de ms 
son bastante similares a los de m, con lo cual 
ambas estimaciones adquieren confiabilidad, ya 
que proceden de enfoques bastante disímiles. 

En las figuras 2, 3 y 4 se muestra el compor- 
tamiento de los registros de las estaciones Villa 
García, La Villita y Juchipila; la primera presenta 
la tendencia ascendente más suave o ligera, y 
la segunda lo contrario. La estación Juchipila 
requirió un análisis especial, ya que su tenden- 
cia detectada es negativa. Con excepción de la 
estación Villa García, todos los registros que 
mostraron tendencia con la prueba de Mann- 
Kendall resultaron no homogéneos con el test 
de Cramer. Como esta prueba se aplicó con dos 
subperiodos de la mitad del registro, en todas 
las estaciones con tendencia, excepto Juchipila, 
se obtuvieron estadísticos negativos para el 
primer subtramo y positivos para el segundo, 
indicando con ello una media menor que la 
general al inicio y mayor en el segundo subpe- 
riodo. Lo anterior indica tendencia ascendente. 

Al comentar con la oficina de Hidrometeo- 
rología de la Dirección Estatal de la Conagua 
en Zacatecas, la tendencia negativa encontrada 
en la estación Juchipila para su registro de tem- 
peratura máxima del mes de mayo, se recibió 
información relativa a un cambio de ubicación 
en el año 1980, el cual modificó sus condiciones 
circundantes. La aplicación de las tres pruebas 
generales y de las 11 específicas al registro de 
la estación Juchipila por subtramos, 1947-1980 
y 1981-2012, indica que el primer lapso es no 
homogéneo y presenta tendencia según el test 
de Kendall y exceso de variabilidad con base en 
el test de Bartlett. En cambio, el segundo lapso 
es homogéneo de acuerdo con la mayoría de 
las pruebas. Por lo anterior, los resultados de la 
estación Juchipila, mostrados en el cuadro 4, no 
serán tomados en cuenta para la formulación de 
las conclusiones. 

Recientemente, algunos autores como Arora 
et al. (2005) y Singh et al. (2008) prefieren dibujar 
las anomalías de la temperatura en lugar de 
sus valores puntuales, como se muestra para la 
estación Monte Escobedo en la figura 5. 
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Cuadro 3. Resultados del test de Mann-Kendall aplicado a las series indicadas y del ajuste de la regresión lineal para detectar tendencia en los registros de temperatura máxima del mes 


Pend. 
ms 


0.0667 


0.0455 


0.0455 


0.1207 


0.0172 


-0.0310 


R** 


m0 


m0 


mx0 


m0 


mx0 


tc 


2.011 


2.000 


1.999 


2.006 


1.998 


4.089 


2.238 


4.003 


10.332 


-3.586 


0.508 


0.275 


0.453 


0.817 


0.227 


-0.409 


Tendencia lineal de la serie X, 


30.176 


31.925 


31.475 


34.609 


35.246 


41.167 


de mayo de las estaciones climatológicas indicadas del estado de Zacatecas, México. 


0.078 


0.041 


0.053 


0.118 


0.016 


-0.043 


T* 


Si 


Si 


Si 


Si 


No 


Si 


Serie Y”, 
Zs 


3.128 


2.218 


3.644 


4.748 


1.429 


-3.409 


Ss 


375 


375 
630 
655 


265 


-617 


T+* 


Si 


No 


Si 


Si 


No 


Si 


Serie X”, 
Zs 


Ss 


2.677 


= 
OQ 
(69) 


1.684 


285 


2/52 


476 


3.238 


447 


1.299 


241 


-2.313 


-419 


T* 


Si 


Si 


Si 


Si 


No 


Si 


Serie X, 
Zs 


3.471 


2.301 


3.262 


6.599 


1.916 


-3.204 


Ss 


416 


389 


564 


910 
1009 
-580 


Nombre 


Monte Escobedo 


Villa García 


San Pedro Piedra G. 


La Villita 


Excamé 


Juchipila 


8 


9 
10 


12 


15 


*Tendencia. **Resultado. 


Respecto al cociente de temperaturas 
máximas medias 


Hu y Liu (2010), en un estudio sobre el cambio 
climático, analizaron 31 registros de lluvia anual 
de amplio periodo (más de 50 años) a través de 
China, encontrado que el cociente del valor 
medio reciente de los últimos 10 años al de lar- 
go periodo, fluctuó de 0.89 a 1.11, con un valor 
medio igual a la unidad. Por ello concluyen que 
la lluvia en tal país no muestra efectos cuantita- 
tivos del cambio climático. 

El cuadro 4 está formado por tres porcio- 
nes. En la primera se exponen los parámetros 
estadísticos de la temperatura máxima del mes 
de mayo: media (O), desviación estándar (S) y 
coeficientes de variación (Cv), asimetría (Cs), 
curtosis (Ck) y correlación serial de orden uno 
(r,). En la segunda parte se tienen las medias (x,) 
de los últimos 5, 10 y 15 años del registro. En la 
tercera se presenan los cocientes de x, entre X. 

Los últimos cinco renglones del cuadro 4 
muestran gran similitud en los indicadores de 
los tres cocientes establecidos, con mayor dis- 
persión para el cociente CO, y una semejanza 
idéntica para los estadísticos de tendencia cen- 
tral de los cocientes CO,, y CO,,. La mediana y 
la moda se estimaron con base en la distribución 
gamma mixta, sin incluir los resultados de la 
estación Juchipila. 

En la columna final del cuadro 4 se tienen 
las diferencias entre la media de un lapso final 
de 15 años y la de todo el registro. Para las 
cuatro estaciones climatológicas que mostraron 
tendencia ascendente —Monte Escobedo, Villa 
García, San Pedro Piedra Gorda y La Villita, 
todas ellas de la Región Hidrológica 12 Parcial 
(Río Santiago)— se tiene una diferencia máxi- 
ma de 2.5 *C y mínima de 1.7 *C, con un valor 
promedio aproximado de 2.1 *C, que puede ser 
considerado el incremento medio que ha mos- 
trado esta subregión del estado de Zacatecas 
debido al cambio climático. 


Impacto del desarrollo urbano 


Seis de las 16 estaciones climatológicas procesa- 
das (cuadro 1) se ubican dentro de áreas urba- 
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Figura 2. Serie cronológica de la temperatura máxima (*C) del mes de mayo de la estación climatológica Villa García, 


Zacatecas, México. 


y = 34.609 + 0.118-x 
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Figura 3. Serie cronológica de la temperatura máxima (?C) del mes de mayo de la estación climatológica La Villita, 
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Figura 4. Serie cronológica de la temperatura máxima (*C) del mes de mayo de la estación climatológica Juchipila, 
Zacatecas, México. 


Pendiente (m) = 0.078 


Anomalía de la temperatura máxima del mes de mayo en “C 


Figura 5. Anomalías de la temperatura máxima ("C) del mes de mayo de la estación climatológica Monte Escobedo, 
Zacatecas, México. 
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nas: Chalchihuites, Zacatecas, Monte Escobedo, 
Villa García, San Pedro Piedra Gorda y Presa 
El Chique. Sin embargo, sólo Zacatecas está 
dentro de una ciudad densamente poblada y el 
resto se localiza en pueblos que incluso tienen 
huertas y áreas de cultivo cercanas, además de 
un escaso número de vehículos automotores, de 
manera que el impacto de la urbanización es 
aún leve o escaso. Resulta interesante destacar 
que en la última columna del cuadro 4, todas 
las estaciones citadas muestran una diferencia 
positiva entre la media reciente y la de amplio 
periodo, destacando San Pedro Piedra Gorda, 
Villa García y Monte Escobedo. 


Conclusiones 


El resultado fundamental de este estudio es que 
la temperatura máxima del mes de mayo, como 
indicadora probable del comportamiento gene- 
ral de la temperatura en el estado de Zacatecas, 
México, está aumentando en la subregión que 
corresponde a la Región Hidrológica 12 Parcial 
(Río Santiago), en la cual cuatro estaciones cli- 
matológicas de las ocho procesadas mostraron 
tendencia ascendente significativa, con un valor 
medio del incremento de 2.1 *C, entre la media 
reciente y la de largo periodo (ver cuadro 4). 

Por lo anterior, se sugiere realizar este aná- 
lisis en otros estados de la república mexicana, 
o bien en otras regiones o subregiones hidro- 
lógicas, para establecer si este aumento de la 
temperatura es generalizado y donde se detecte, 
como en la subregión citada del estado de Zaca- 
tecas, México, establecer medidas orientadas a 
redimensionar los sistemas de abastecimiento 
de todo tipo, pues la evaporación y la evapo- 
transpiración sin duda estarán aumentando, 
disminuyendo con ello la disponibilidad en 
los embalses e incrementando las necesidades 
hídricas de los cultivos tanto bajo riego como 
de temporal. 
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Apéndice 1: Referencias para 13 de las 16 
pruebas estadísticas aplicadas 


1. Test del Cociente de Geary (CG) 


Relaciona la desviación media absoluta del 
registro con su desviación estándar. Cuando 
el CG se aparta de la unidad, indica pérdida 
de normalidad. Su expresión es la siguiente 
(Machiwal €z Jha, 2012): 


(A.1) 


2. Test de Shapiro-Wilk (SW) 


También conocido como Test W, emplea valores 
tabulados para calcular su estadístico. Su proce- 
dimiento y datos necesarios se pueden consultar 
en Shapiro y Wilk (1965), y en Ruiz-Maya (1977) 
para registros de hasta 50 datos, y en Shapiro 
(1998) para series de hasta 98 valores. 


3. Prueba de Helmert (PH) 


Esta es una prueba no paramétrica que se basa 
en secuencias y cambios respecto a la media, 
la cual se pueden consultar en Mather (1977) y 
Campos-Aranda (1998). 


4. Prueba de las secuencias (PS) 

Esta prueba no paramétrica emplea las secuen- 
cias que se forman respecto a la mediana; se 
pueden consultar en Thom (1971) y Campos- 
Aranda (1998). 


5. Prueba de Von Neumann (PV) 


Esta prueba paramétrica detecta pérdida de 
aleatoriedad por componentes determinísticas 
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no especificadas; su procedimiento operativo 
se puede consultar en WMO (1971), Buishand 
(1982), y Machiwal y Jha (2012). 


6. Test de Anderson de persistencia (TA) 


Esta es una prueba paramétrica que se basa en 
el coeficiente de correlación serial de orden 1 
(r,), la cual se pueden consultar en WMO (1971), 
Linsley, Kohler y Paulhus (1988), y Campos- 
Aranda (1998). 


7. Test de Sneyers de persistencia (TS) 


Es una prueba no paramétrica, recomendada 
cuando el registro no es normal. Se puede con- 
sultar en WMO (1971), y Machiwal y Jha (2012). 


11. Test de Bartlett de variabilidad (TB) 


Esta prueba paramétrica, también llamada de 
homogeneidad de varianzas, busca inconsistencia 
debido a la dispersión. Su secuela se puede 
consultar en WMO (1971) y Ruiz-Maya (1977). 


12. Test de Cramer de cambio en la media (TC) 


Esta prueba no paramétrica compara la media 
de un subperiodo con la del registro; se puede 
consultar en WMO (1971) y Campos-Aranda 
(1998). 


13-16. Tests de Buishand de cambio en la 
media (TB, a TB,) 


La gráfica de las desviaciones acumuladas 
con respecto de la media se denomina curva 
masa residual, permite reconocer fácilmente 
los cambios que han ocurrido en ésta y por 
ello Buishand (1982) presentó cuatro pruebas 
basadas en estadísticos (Q, R, U y A) tomados 
de tal curva. Aplicaciones de estas pruebas en 
México se pueden consultar en Gutiérrez-López, 
Arellano-Lara y Aparicio-Mijares (2007), y en 
Campos-Aranda (2008). 


Apéndice 2: Descripción de los tres test 
de tendencia 


Machiwal y Jha (2008) aplicaron 12 pruebas para 
detección de la tendencia, la mayoría paramé- 
tricas y encontraron que los tres test siguientes 
fueron los de mayor potencia o más confiables. 


8. Test de Kendall (TK) 


Esta prueba paramétrica también se conoce 
como de correlación de rangos (Machiwal « Jha, 
2012). Es muy efectiva si la tendencia funda- 
mental es de tipo lineal o se aproxima a ésta. Su 
procedimiento consiste en evaluar el número p 
en todos los pares de observaciones (x, x; con 
¡=1lan-1;j=i+1an) en que x, es mayor que 
x; después se calculan las estadísticas siguientes 
(Kottegoda, 1980): 


7 j 
ee (A.2) 
, 2(2n+5) 

5 9n(n-1) o 


Si el cociente 1/o, es menor que +1.960, el 
registro no presenta tendencia en una prueba de 
dos colas con un nivel de significancia a de 5%. 


9. Test de Spearman (TS,) 


En esta prueba no paramétrica, que se conoce 
como correlación del orden de rangos (Machiwal 
é£z Jha, 2012), primero se ordenan los datos de 
menor a mayor y se sustituye cada elemento 
del registro por su rango (k,); después, para 
cada dato de la serie original, se calculan las 
diferencias: d,=k, - 1, coni=1a n y se evalúan 
las estadísticas siguientes (WMO, 1971; Adeloye 
éz Montaseri, 2002): 


6 Ya 
Risle= A 


an (A.4) 
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(A.5) 


El valor de la ecuación anterior se compara 
con el intervalo que define +f_ tomado de la 
distribución t de Student con n - 2 grados de 
libertad y nivel de significancia a, para una 
prueba de dos colas; si lo excede (t > t.), existe 
tendencia. 


10. Test de Mann-Kendall (TM) 


Es una prueba no paramétrica que busca ten- 
dencia en la serie sin importar si es lineal o no 
lineal; pero el registro no debe mostrar persis- 
tencia. Su estadístico operativo Ss es la suma 
de los signos de las diferencias (z) de todas las 
parejas factibles de formar, por ello su expresión 
es (Manly, 2001): 


n il 


Ss = Y Y signo(x, -x,) 


122 jul 


(A.6) 


en la cual el signo (z) es -1 para z< 0, O para 
z=0 y +1 para z> 0. Para una serie con valores 
aleatorios se espera que Ss = 0, cuya varianza es 
(Hirsch, Helsel, Cohn, éz Gilroy, 1993; Machiwal 
é£z Jha, 2012): 


n(n -1)(2n+5) 


Var(Ss) = (A.7) 
El estadístico de la prueba es: 
Z¿= 54m, (A.8) 
ar(Ss) 


En la expresión anterior, m_=1 cuando Ss <0 
y m,=-1 para Ss > 0. Si el valor absoluto de Z, 
es mayor que el valor crítico de la distribución 
Normal estándar, la serie presenta tendencia 
creciente o decreciente con un nivel de signifi- 
cancia a. Para a, = 5%, el valor crítico es 1.960. 
La varianza (ecuación (A.7)) se debe corregir 
cuando existen datos iguales o la muestra es 
censurada (Hirsch ef al., 1993). 


Hamed y Rao (1998) desarrollaron una 
corrección para la ecuación (A.7), que toma en 
cuenta la persistencia de la serie, la cual permite 
la aplicación del Test de Mann-Kendall a datos 
autocorrelacionados. Sonali y Nagesh-Kumar 
(2013) presentan dos enfoques de corrección de 
la prueba de Mann-Kendall basados en modi- 
ficar su varianza (ecuación (A.7)). Kundzewicz 
y Robson (2004) exponen otros procedimientos 
de aplicación de la prueba de Mann-Kendall, los 
cuales utilizan el re-muestreo estadístico. 


Apéndice 3: eliminación de la 
persistencia por pre-blanqueado del 
registro 


Desde mediados de la década de 1990 se 
estableció que la correlación serial positiva 
de una muestra o registro incrementa la posi- 
bilidad de que el Test de Mann-Kendall (TM) 
detecte tendencia significativa. Para eliminar 
tal influencia se propuso el pre-blanqueado del 
registro, considerando que tal componente de 
correlación serial es un proceso autorregresivo 
de orden uno [AR(1)] y que por lo tanto debe 
ser retirado antes de aplicar el test TM (Burn éz 
Hag-Elnur, 2002). 

Posteriormente, Yue, Pilon, Phinney y 
Cavadias (2002) demostraron, mediante si- 
mulación aleatoria, que cuando un proceso 
AR(Í) y una tendencia lineal existen en una 
serie, la remoción del proceso autorregresivo 
con el pre-blanqueado también elimina una 
parte de la tendencia. Lo anterior se observa 
considerando que un registro (X,) está formado 
por una tendencia lineal (T, = f$*t) y un proceso 
AR(1), con r, como el coeficiente de correlación 
serial de orden uno y e, como el ruido blanco o 
componente aleatoria (A, =1¡A,+€ clk es decir: 


LEA (A.9) 
El registro preblanqueado será: 
X, =X, 5 X,, (A.10) 


Sustituyendo en la expresión anterior la 
ecuación (A.9) y el proceso AR(1) se tiene: 
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X, =B-*+A,-1 [B(t-1)+4,,] 


=p-t+(r, A, +6) 1 Plt-1)-7,:4,, (AM) 
Simplificando y agrupando se tiene: 
X, =(1-1,):B:t+8,+1,:P (A.12) 


En la expresión anterior se observa que el 
proceso autorregresivo ha sido eliminado, pero 
ahora la pendiente de la serie pre-blanqueada es 
P' = (1- r,) : P; entonces, cuando se remueve un 
AR(1) positivo (r, > 0) se tiene que If"! < IfBl, 
indicando que el pre-blanqueado elimina parte 
de tendencia. Cuando AR(1) es negativo (r, < 0), 
ocurre lo contrario, |1f'l > 1f 1, de manera que 
el pre-blanqueado incrementa la tendencia. 

Debido a lo anterior, Yue et al. (2002), y 
Yue, Pilon y Phinney (2003) propusieron un 
pre-blanqueado libre de tendencia, que consiste 
en identificar la tendencia lineal (T,) y retirarla 
de la serie (X,), para obtener una muestra sin 
tendencia (Y), esto es: 

Y =X, -T,=X,- Pt (A.13) 

Después se elimina el proceso AR(1), cuando 
el r, resultó estadísticamente significativo o di- 
ferente de cero, para obtener la serie de valores 
independientes (Y), ésta es: 

de La (A.14) 

Finalmente, se suma la tendencia identifi- 
cada para obtener una muestra mezclada de 
que incluye tendencia y otros componentes, 
como el ruido blanco, pero que no está influen- 
ciada por la correlación serial; a esta serie se le 
aplica el test de Mann-Kendall: 

YY AT Y + po (4.15) 
Apéndice 4: Recta de regresión para la 
tendencia lineal 


Se considera que la variable dependiente (y) son 
las temperaturas anuales T, en *C y los tiempos 


o años A, son las abscisas (x), en este caso iguales 
al ¡-ésimo valor ¡. Para probar si la pendiente (m) 
de la recta de regresión ajustada por mínimos 
cuadrados de los residuos es estadísticamente 
diferente de cero, se usa una prueba basada 
en la t de Student definida por las ecuaciones 
siguientes (Ostle $ Mensing, 1975): 


T=b+m:A (A.16) 
b=T-m-:A (A.17) 
1 T.A 
Cov(T,A) 2 dd 
7 Var(A) = E a (A.18) 
yA -A 
0 Cov(T,A) (4.19) 
Y NVar(T)-Var(A) 
Var(T)= ST ed (A.20) 
i=1 
is (4.21) 
S 
2: 
=> 5 : (A.22) 
$41 
i=1 
(12) 
2 — tel 
S? = 72) (A.23) 


T y A son las medias aritméticas; r,, es el 
coeficiente de correlación, el cual varía de cero 
a +1. La pendiente mm tiene unidades de *C/año 
y la ordenada al origen b de *C. En la ecuación 
(A.23), T es el valor estimado con la recta de re- 
gresión (ecuación (A.16)). S? y S? son las varian- 
zas de los errores y de la pendiente de la recta 
de regresión. Si el valor calculado t (ecuación 
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(A.21)) es mayor que el crítico (t ), obtenido para 
la distribución t de Student con n - 2 grados de 
libertad y a = 5%, en una prueba de dos colas, 
la pendiente m es significativa, es decir, existe 
tendencia lineal. El problema de esta prueba es 
que no distingue entre persistencia y tendencia 
(Adeloye €: Montaseri, 2002); además, la estima- 
ción de m es sensible a la no normalidad de los 
datos (Sen, 1968). 


Apéndice 5: Pendiente de la tendencia 
(ms) según criterio de Pranab Kumar Sen 


Sen (1968) propuso utilizar el criterio de análisis 
de la tendencia establecido en el test de Kendall 
para estimar su pendiente (ms), definiéndola 
como el valor mediano (MED) de las pendientes 
parciales, es decir: 


XxX -X del 
es para toda j¡>i (A.24) 
j-i 


X, y X, son los datos, en este caso temperaturas 
máximas de mayo, en los tiempos j¡ e í, que son 
años; ms tiene unidades de *C/año. Siendo n 
el número de datos, entonces habrá nín - 1)/2 
pendientes parciales. El signo positivo de ms 
define tendencias ascendentes y el negativo 
descendentes (Sonali £ Nagesh-Kumar, 2013). 
Machiwal y Jha (2012) indican que la estima- 
ción ms es resistente o robusta a la presencia de 
valores dispersos (outliers). 
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